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FOREWORD 


The  work  described  in  this  report  was  performed  in  the  Fire  Control  Formulation  Branch 
(K41),  Submarine  Launched  Ballistic  Missile  (SLBM)  Research  and  Analysis  Division,  Strategic 
and  Space  Systems  Department,  and  was  authorized  under  Strategic  Systems  Program  Office  Task 
Assignment  36401.  This  work  was  necessitated  by  the  need  to  formulate  an  exact  method  of 
computing  geodetic  latitude  and  altitude  whose  chief  attribute  is  the  absence  of  singularities  in  the 
results. 


This  report  is  a  refined  version  of  NSWC  TR  85-85,  dated  April  1985  and  revised  May 
1986.  NSWC  TR  85-85  should  be  discarded. 

This  report  has  been  reviewed  and  approved  by  Davis  Owen;  Johnny  Boyles;  Carol  Rose, 
Head,  Fire  Control  Formulation  Branch;  and  Dr.  David  Lando,  Head,  SLBM  Research  and 
Analysis  Division. 


Approved  by: 


R.  L.  SCHMIDT,  Head 

Strategic  and  Space  Systems  Department 
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INTRODUCTION 


A  simple  and  efficient  model  is  presented  for  calculating  the  exact  geodetic  latitude  and  altitude 
of  an  arbitrary  point  in  space,  given  the  coordinates  of  that  point.  The  model  assumes  the  Earth  to 
be  an  oblate  spheroid;  i.e.,  an  ellipsoid  of  revolution  whose  semimajor  axis  ‘a’  is  the  radius  of  the 
circle  described  by  the  equatorial  plane,  and  whose  semiminor  axis  ‘b’  is  a  line  joining  its  center 
and  one  of  its  poles.  The  point  in  question  can  be  either  inside  or  outside  the  Earth  ellipsoid, 
excluding  a  small  open  region  bounded  by  a  prolate  spheroid  that  is  concentric  with  and  contained 
within  the  Earth  ellipsoid.  The  cardinal  merit  of  this  model  is  the  absence  of  singularities  in  the 
results. 


FORMULATION 


Choose  an  Earth-fixed  Cartesian  coordinate  system  whose  origin  coincides  with  the  center  of 
the  Earth  ellipsoid.  The  unit  vectors  i,  j ,  and  Ic  coincide  with  the  x,  y,  and  z  axes,  respectively. 
The  +z  axis  points  in  the  direction  of  the  North  Pole.  The  +x  axis  is  the  line  of  intersection  of  the 
equatorial  plane  with  the  plane  of  zero  longitude.  The  +y  axis  completes  a  right-handed  coordinate 
system. 

An  equation  of  the  ellipsoid  in  this  frame  is 

-H  -^  =  1,  where  r  =  Jx^  +  y^  . 
sr  b^ 

Let  P(xo,  yo,  zq)  be  the  coordinates  of  the  given  point.  It  is  desired  to  find  the  nearest  point 
PXx,y,z)  to  P  on  the  surface  of  the  ellipsoid. 

The  slope  of  a  normal  to  the  ellipsoid  at  any  point  on  its  surface  is  given  by 

1  a^z  ,  r^  z^ 

--7-  =  where  -5- --7  =  1 .  (2) 

^  b^r  a^  b^  ^  ^ 

dr 

Therefore,  the  slope  of  a  normal  from  P  is 

— - =  rr^  where  ro  =  ^Xq  -t-  yo  ; 

ro-r  b"r 

i.e., 

(ro-r)a'z  =  (zo-z)b'r. 
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or 


r.a^z  =  r{(a'-b^)z  +  b%}. 

Squaring  both  sides  and  expressing  r  in  terms  of  z, 

aVr,V=(b^-z^){(a^-b=)z+b^z.f. 

Writing  the  above  equation  in  descending  powers  of  z,  the  following  quartic  in  z  is  obtained: 

(a’  -  b^f  z^  +  2b“(a^  -  b^)z.z’  +  b^{aV^  +  b^zj  -  (a^  -  b^f  }z^  -  2b‘(a^  -  b^)z.z-  b^zj  =  0  . 

Since  z  is  either  positive  or  negative  and  sign  z  =  sign  Zq,  the  above  quartic  can  be  expressed  in 
terms  of  z  .  The  result  is 


zf  +  2p|zp  4-  q|zf  -  2pb^|z|  -  p^b^  =  0  , 


(3) 


where 


b^z, 


0  _  ro 


P  a^-b"  e'^  ’ 


(4) 


b^{a\Hb^zS-(a^-b^f} 

q  =  - 


^2  .  2  . 
p  -  b  + 


»2„/2  ’ 
e  e 


(5) 


where 


e"  =  l-^,e'^=4-l 

a  b^ 


Two  of  three  powerful  theorems  in  the  Theory  of  Equations  will  now  be  employed  to  expose 
the  nature  of  the  roots  of  Equation  (3).  The  third  theorem  will  be  needed  later  on.  These  theorems 
are  found  in  the  referencei. 


I.  An  equation  f(x)  =  0  cannot  have  more  positive  roots  than  there  are  changes  of  sign  in  f(x), 
and  cannot  have  more  negative  roots  than  there  are  changes  of  sign  in  f(-x). 

n.  Every  ^nation  that  is  of  an  even  degree  and  has  its  last  term  negative  has  at  least  two  real 
roots,  one  positive  and  one  negative. 

in.  Every  equation  of  an  odd  degree  has  at  least  one  real  root  whose  sign  is  opposite  to  that  of 
its  last  term. 


By  observing  Equations  (3),  (4),  and  (5),  the  following  conclusions  are  deduced: 


1  Hall,  H.S.  and  Knight,  S.R.,  Higher  Algebra.  Fourth  Edition,  MacMillan  and  Co., 
London,  1927. 
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(i)  Since  the  last  term  of  Equation  (3)  is  negative,  there  are  at  least  two  real  roots  of  opposite 
sign  (using  II). 

(ii)  Since  there  is  only  one  change  of  sign  in  Equation  (3),  there  is  at  most  one  positive  root 
(using  I). 

From  (i)  and  (ii),  it  is  seen  immediately  that  Equation  (3)  has  exactly  one  positive  root,  and  this 
is  the  root  that  is  sought. 

The  solution  of  Equation  (3)  is  effected  by  a  standard  method  known  as  Ferrari’s  method^ .  In 
the  equation 


|zf  +  2p|zp  +  q|z|^  -  2pb^|z|  -  p^b^  =  0  , 

add  to  each  side  (clzl+d)2,  the  quantities  c  and  d  being  determined  so  as  to  make  the  left-hand  side 
a  perfect  square;  then 

\zf  +  2plzf  +  (q  -I-  c^)|zf  -I-  2(cd  -  pb^)jz|  +  d^-  p^b^  =  (c|z|  +  d)^  . 


Suppose  that  the  left  side  of  the  equation  is  equal  to  (Izl^  +  plzl  + 1)^,  then  comparing  the 
coefficients. 


p^  -h  2t  =  q  +  c^,  pt  =  cd  -  pb^  t^  =  d^  -  p^b^ . 
Eliminating  c  and  d  from  these  equations, 

p^(t  +  b“f=(2t  +  p^-q)(t^+pV). 


(6) 


or 


or 


2t^  -  qt^  -  p V(q  -  p^  -I-  b^)  =  0  , 


2t^  -  qt^  -  ^  ^  =  0,  using  Equation  (5) 

e  e 


or 


2^2„2 


2e-qe-^^=o 


(7) 


Applying  Theorems  I  and  m  to  Equation  (7),  we  see  that  Equation  (7)  has  exactly  one  positive 
root.  The  solution  of  Equation  (7)  is  accomplished  by  a  standard  method  known  as  Cardan’s 
Solution!.  Eliminating  the  t2  term  in  Equation  (7)  by  making  the  substitution 
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f  nV  f  nV 

2t'+s  -q,'+a  -ai^=o, 

\  6J  y  6J  e 


where 


t'^+Ft'  +  H  =  0  , 


JP  _  _  T.  pj  _  _ ^  ^0  ^0 

“  12’  “  108  2e'* 


The  solution  of  Equation  (8)  is 


,  H 

^  “  n  4  ■‘'27  2 


F^  H 
4  27  2 


Equation  (9)  is  valid  provided  that 


We  have 


where 


4  27 


aV7-2  9 

H  =  — 3 - =  _  2  /  3^2P), 

108  2e'*  6^^  ’ 


108  2e' 


27^ 

e'* 


F^  1 


^^>0. 

4  27 


P  +  q'>0, 

which  is  satisfied  when  q  >  0;  i.e.,  when  a^r^  +  ^  (a^  -  )  ,  which  represents  all  space 

excluding  the  open  region  bounded  by  the  ellipsoid 
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2  2 
aV  bV" 


Since  ae2  <  be '2  <  43  km,  this  region  is  of  no  practical  interest.  Therefore,  q  ^  0  will  be  the 
constraint.  Applying  Theorem  I  to  ^nation  (7)  with  -t  substituted  for  t,  we  see  that  Equation  (7) 
has  no  negative  roots.  Hence,  Equation  (7)  has  only  one  real  root  and  it  is  positive. 


I  4  27  2  “  6' 


where 


Similarly, 


where 


/P(P+ q’)  +  ^(2P+ q’)  =  ^[2P+ q=  +  2^P(P  +  q=)} 


f /h^ 

.Hi 

a 

^  27 

2  ) 

Q=(VP+q’  +Vp)'  ■ 


U  27  2  ,  “  6 


q'  =  (Vp+?-Vp)'  . 


(Note  that  ^  “  q  *  However,  putting  Q'  in  this  form  is  undesirable  as  it  introduces  a  singularity 
when  Q  =  0,  which  occurs  when  either  both  q  and  ro  equal  zero  or  both  q  and  Zq  equal  zero.) 
Substituting  these  values  into  Equation  (9), 

t'  =  7(Q+Q')=>i  =  7(q+Q+Q')- 

O  O 

It  will  now  be  shown  that  q  <  2t  <  q  +  b2,  where  the  left-hand  equality  occurs  when  either 
ro  =  0  or  Zo  =  0  and  the  right-hand  equality  occurs  when  aro  =  bizol. 


Rewriting  Equation  (7)  in  the  form 


^2^2  2 
3.  Tn  Zn 


t^(2t-q)  =  ^, 

e 


the  lower  bound  of  t  is  deduced  at  once.  Now  let 
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where  t  >  q.  Then 


3  2 

g  =  T^-qx^ - ^ 

e 


=  t(3x  -  2q)  >  0  . 


Thus,  g  is  a  monotonically  increasing  function  of  x  for  all  x  >  q.  Since  2t  >  q  and  g  =  0  when 
X  =  2t,  it  follows  that  g  >  0  when  x  >  2t.  Putting  x  =  q  +  in  the  expression  for  g. 


4d  p  n 

— ^sing  Equation  (5) 
e  e 


I  e  e  J 


Tq  Tq  Zq 

2t<q  +  b^,  equality  occurring  when  p= — i.e.,  when —  =  — r,  or  when  arQ  =  b|zg|,  in 

CC  C  C  '  ' 


which  case, 


q  =  2p^-b^  t  =  |(q  +  b^)  =  p^ . 


Solving  for  c  and  d  from  the  system  of  equations  (6), 

c  =  -7p^-q  +  2t,  d  =  -(t  +  b^)  =  -^/t^  +  pV  . 
c 

Since  2t  >  q,  it  follows  that  c  >  p.  Now,  (Iz|2  +  plzl  + 1)2  =  (clzl  +  d)2  implies 
Izl2  +  plzl  + 1  =  ±  (clzl  +  d),  from  which  the  following  two  quadratics  in  Izl  are  obtained: 

jzf +  (p-c)|z|  +  t-d  =  0  , 


|zp+(p  +  c)|z|  +  t  +  d  =  0 .  (12) 

It  has  already  been  shown  that  Equation  (3)  has  exactly  one  positive  root.  Furthermore,  since 
t2  -  d2  =  -p2b2  <  0  and  d  >  0,  it  follows  that  t-d<0<t  +  d.  Now,  t  -  d  is  the  last  term  of 
Equation  (11)  and  has  just  been  shown  to  be  negative.  Hence,  by  applying  Theorem  n  to  Equation 
(1 1),  we  see  immediately  that  Equation  (11)  must  contain  the  required  positive  root.  Hence, 
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-(p-c)  +  v(p-cf -4(t-d) 
2 


=  — l^c  -  p  +  ^j2p^  -  q  -  2t  -  2pc  +  4^ +  p^b^  J  . 

This  solution  is  free  of  singularities.  It  is  easy  to  show  that  lim|z|  =  0  and  lim|z|  =  b  .  We  have 

Zo->0  Iq^O 

lim  p  =  0,  lim  P  =  lim  P  =  0  =>  lim  Q  =  lim  Q'  =  lim  Q  =  lim  Q'  =  q 

Zo->0  Zq-^O  Zo->0  Zo-»0  Tq-^O  Iq-^O 


=>  lim  t  =  lim  t  =  —  lim  c  =  0,  lim  c  =  p. 

^0-40  2  '0*^0 


To  compute  the  geodetic  latitude  and  altitude  at  P'(x,y,z),  it  is  desirable  to  introduce  a 
geometric  term,  Ne,  which  is  never  zero.  Ng  is  defined  to  be  the  distance  along  the  ellipsoidal 
normal  from  the  surface  of  the  ellipsoid  to  the  z-axis  (see  Figure  1). 
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z 


HGUREl.  ELLIPSOIDAL  NORMAL 


From  Figure  1, 


Also, 


3.  Z 

land)  =  —5-,  from  Equation  (2) 
b  r 


=rsec<j)  =  r-(/l  +  tai^  =  + 


and 


,4,2 


z-  1  a'z'  L  .  a  -b  2 


b^ 


—  1  +  ' 


=aj] 


a^z 

sind)  =  cos<l)tan(t)  =  — r — 
b  N. 


•  (j)  =  sin 


■1  a^z 
b^N. 


which  is  the  desired  expression  for  the  geodetic  latitude. 
Also  from  Figure  1, 


h  cost])  =  To  -  r  =  ro  -  cos(t). 


NSWCDD/rR-94/77 


hsin<t)  =  Zq  -  z  =  Zo  — r-N^sin<l) 


Multiplying  the  first  equation  by  cos  ^  and  the  second  by  sin  ^  and  adding, 


h  =  roCos(l>  +  ZoSin(|)-N^  cos  ([)  +— sin  <])  . 


.2  a  , 


cos^<t)  +  -5-sin^<l): 


2  .  a'^z^ 

tan  <|)  =  ^4^2 

a^z^ 

■b^(b^-z^) 

1  b"  2a 
a 

^  ^  u2  2 

_  D  — Z 

1  +  tan^  ([) 

2  2 
a  z 

1  *•  *>  /  »!  ' 

b^(b^-z^) 


b‘+(a^-b^)z= 


Hence 


1  a 

,2  2  =  — 2  »  using  Equation  (13). 
e  z  N: 


h  =  Tq  cos<t>  +  Zq  sin  ((>  -  —  , 


which  is  the  desired  expression  for  the  geodetic  altitude. 

To  recapitulate,  given  a,  b,  xq,  yo,  and  zq,  the  algorithm  is  as  follows 

7  .  b  ,0  B.  -  I  T~  7  l^ol  7 


p2_,  ^  _ ,  _  /„2  .  2  __N  n-n^  I 

C  —  1  2’®  t.2  ^’^0  •V^O^yO’P  /2  ’  9  P  ®  ^  »2  /2  ‘ 

3.  D  C  C  6 


If  q  >  0 ,  then 


2  2 

p  =  — Q  =  (VP  +  q'+Vp)l  Q'  =  (Vp  +  q'  - a/pP,  t  =  7(q  +  Q  +  Q'), 


t,  z==— |^c-p  +  *y2p^  -  q-2t-2pc  +  4-^/t^H-pV  Jsign  Zq  , 


e'^z^  __  ^2^ 

N^  =  a..  1  +  —J-,  <|)  =  sin  '-n— ,  h  =  roCOS(l)  +  ZoSin(j)-—  . 
V  b  b  N, 
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CONCLUSION 


An  efficient  model  for  transforming  Earth-centered,  Earth-fixed  coordinates  to  geodetic 
coordinates  has  been  presented.  The  model  systematically  derives  exact  expressions  for  the 
geodetic  latitude  and  altitude,  which  are  free  of  singularities.  To  the  author’s  best  knowledge,  no 
such  expressions  exist  in  the  literature. 
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